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Abundance stratification in Type la Supernovae - 
II: The rapidly declining, spectroscopically normal 
SN 2004eo 
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ABSTRACT 

The variation of properties of Type la supernovae, the thermonuclear explosions 
of Chandrasekhar-mass carbon-oxygen white dwarfs, is caused by different nucleosyn- 
thetic outcomes of these explosions, which can be traced from the distribution of abun- 
dances in the ejecta. The composition stratification of the spectroscopically normal 
but rapidly declining SN 2004eo is studied performing spectrum synthesis of a time- 
series of spectra obtained before and after maximum, and of one nebular spectrum 
obtained about eight months later. Early-time spectra indicate that the outer ejecta 
are dominated by oxygen and silicon, and contain other intermediate-mass elements 
(IME), implying that the outer part of the star was subject only to partial burning. 
In the inner part, nuclear statistical equilibrium (NSE) material dominates, but the 
production of 56 Ni was limited to ~ 0.43 ± 0.05Af Q . An innermost zone containing 
~ 0.25Af Q of stable Fe-group material is also present. The relatively small amount 
of NSE material synthesised by SN 2004eo explains both the dimness and the rapidly 
evolving light curve of this SN. 
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1 INTRODUCTION 

Type la supernovae (SNela), the thermonuclear ex- 
plosions of carbon-oxygen white dwarfs that approach 
the Chandrasekhar mass limit (for a review, see 
iHillebrandt fc Niemever]|2000l ). are characterised by a sim- 
ple optical light curve, showing a rise to maximum followed 
by a decline. The light curve is governed by the deposition 
of the 7-rays and the positrons produced by the decay of 
56 Ni, which is copiously produ ced in the explosion and de- 
cays into 56 Co and hence 56 Fe (|Arnettlll982l ; lKuchner et all 
1994), and by the diffusion in the expanding ejecta of the 
optical photons that are produced in the thermalisation pro- 
cess. 

A correlation between t he peak luminosity of SNe la and 
the shape of the light curve (|Phillipslll993l : lRiess et al.lll996l : 
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iPhillips et allll999l h makes this class of SNe good standard- 
izable candles. SNela have been used to explore the expan- 
sion propert ies of the Universe out to about one third of its 
present age l|Riess et al.| [T998: Pcrl mutter et al.1l 19971 ). 

While the SN luminosi ty is related fairly directly to 
the amount of 56 Ni ejected (|Arnettl 1 19*821) ■ so that the ob- 
served spread in peak luminosity suggests a correspond- 
ing spread in 56 Ni mass of about one order of magnitude 
f~ 0.1 to ~ lMnl fContardo. Leibundgut, fc Vaccal |2000| ; 



IStritzinger et alj|200tt iMazzah et al.ll2007ah the light curve 
width depends primarily on the optical opacity in the ejecta, 
which is a function of both temperature and composition 
jHoeflich et alJll993h . In the H-free ejecta of SNela, line 
opaci t y dominates over continuum processes dKarp et al.l 
1 19771 : iHoeflich et al. I ll993l : iPauldrach et al.lll996l ). Heavier 
elements have a thicker array of lines which, in combina- 
tion with the rapid and differential expansion of the SN 
ejecta, causes the effective frequencies of the lines to overlap 
over the volume of the ejecta, efficiently blocking most UV 
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and blue photons. For most of these photons the only way 
to escape is to be re-emitted in some less saturated tran - 
sition in the red l|Pinto fe Eastman! l2000l : iMazzalil l2000h . 
This "fluorescence" process determines the shape of Type 
I SN s pectra at early times , when the density is sufficiently 
large. iMazzali et al.l (|200ll ) showed that if the opacity is 
parametrised according to the average number of line tran- 
sitions of the species that dominate the ejecta (singly and 
doubly ionised Si, S, Fe, Co and Ni for SNefa), SNela that 
produce more 56 Ni are not only brighter but also have a 
larger opacity and hence broader light curves. 

Additionally, different conditions within the white 
dwarf (e.g. metallicity) c an lead to different nucleosynthe - 
sis of Fe- group isotopes l|Timmes. Brown, fe Truranl 12003). 
For the same production of NSE material, SNe with 
a higher fraction of stable isotopes (e.g. 54 Fe, 58 Ni) 
over radioactive 56 Ni have light curves with a simi- 
lar shape but different brightness. This may cause a 
spread in the observed luminosi ty-light curve shape relation 
(|Mazzali fe Podsiadlowskll2006l ). 

Although an explanation for the observed behaviour of 
SNela seems to have been found, a conclusive proof of the 
nature of the physical mechanism behind these explosions 
still eludes us. Subsonic deflagrations of Chandrasekhar- 
mass carbon-oxygen white dwarfs seem to be unable to 
prod uce the required energy and 56 Ni yield l|Ropke et al.l 
120071 ). Detonations incinerate the white dwarf to NSE, and 
do not produce the IME that are observed in t he spectra. 
A tra nsition from a deflagrat ion to a detonation |Khokhlovl 
Il99ll ) is a possible solution l|Mazzali et al.l l2007ah . but its 
physical basis remains to be established, despite much effort 
jRopke fe Niemeverll2007l ; lR5pkell2007l : IWooslevll2007h . 

An alternative approach to the study of SNe la, rather 
than computing explosions and predicting their appearance, 
is to map the outcome of real events. This can be done per- 
forming spectrum synthesis of a time-series of SN spectra, 
following the progressive exposure of deeper and deeper lay- 
ers of the ejecta as they expand. 

Spectra at early times, in the so-called photospheric 
phase, probe only the outer half or so of the depth of the 
ejecta. A super-nebular phase then follows, during which 
the position of the pseudo-photosphere is extremely wave- 
length dependent and some features already follow nebu- 
lar physics. A more direct and complete view of the inner 
ejecta can be obtained about one year after the explosion, 
when the gas in the young SN remnant is heated by colli- 
sions with the fast particles from the deposition of 7-rays 
and positrons and cools by emission of radiation, mostly in 
forbidden [Fe 11] and [Fe ill] lines. In this phase, the domi- 
nance of Fe lines testifies that 56 Ni is synthesised mostly in 
the deepest p arts of the star , in lin e with theoretical expec- 
tations (e.g. Ilwamoto et al.l Tl999) . A relation between the 
width of the Fe emission lines - which indicates the velocity 
of expansion of the bulk of 56 Ni - and both the SN peak 
brigh tness and the light curve decline rate l|Mazzali et al.l 
119981 ) confirms that in more luminous SNe a larger fraction 
of the star is burned to NSE. It is therefore essential that 
SNe that have been monitored in the peak phase are also ob- 
served spectroscopically and photometrically in the nebular 
phase. 

Only few examples of S Nela with very de tailed time 
coverage were available (e.g. ISalvo et al.ll200lh before the 



European Research and Training Network on the physics of 
Type la Supernovae (RTN) began collecting dat a of nearby 
SNela systematically (e.g. iBenetti et al.l 12004 ). Now, the 
availability of a sizable sample of SNela with detailed ob- 
servational coverage makes it possible to study in greater 
depth the properties of SNe of different luminosities. Since 
such differences are the result of different nucleosynthetic 
outcomes of the explosions, deriving the composition of the 
ejecta from the data should shed light on the details of the 
explosion process and it variations. 

A fter initial efforts u sing homogeneous compositions 
(e.g. IMazzali et ail Il993l ). the more computationally ex- 
pensive analysis of a time series of spectra using stratified 
abundances has b een attempted in only one case so far. 
IStehle et alj (120051 . hereafter Paper I) mapped the ejecta 
of the normally luminous SN 2002bo from the centre out to 
v ~ 20, 000 kms -1 . They showed that abundances are strat- 
ified but some mixing does occur. They detected burning 
products, namely stable Fe, Si and Ca, at very high veloci- 
ties. This material gives rise to high velocity features (HVF), 
an appa rently ubiquitous pr operty of SNela at the earliest 
phases (jMazzali et ail 120051 ) . For SN2002bo, the IME re- 
gion extends down to ~ 10,000 tons" , below which NSE 
material dominates. 

This distribution of elements in SN 2002bo seemed to 
favour a Delayed Detonation scenario. When the abundances 
derived in this sort of tomographic scan of the ejecta were 
used to generate opacities in a ~ 10 51 erg, Chandrasekhar- 
mass e xplosion with the den sity structure of the classic W7 
model (|Nomoto et ailll984l ). the observed bolometric light 
curve was reproduced extremely well, including the early 
rise. This is induced by mixing out of 56 Ni, which is not 
predicted by most current expl osion models but i s indeed 
detected in early SNIa spectra (|Tanaka et ai|[2008l ). 

In this series of papers our aim is to describe the com- 
position of the ejecta of SNela with different luminosities 
and light curve decline rates. SN 2002bo was a SN la of av- 
erage luminosity, which synthesised ~ 0.5M@ of 56 Ni. In 
this second paper of the series, we address a SN at the dim 
end of the distribution of spectrosc opically normal SNela . 
SN 2004eo was observed by the RTN dPastorello et al.ll2007l ). 
It was a relatively underluminous event, with a B-band post- 
maximum decline rate Ami5(_B)= 1.46 mag, which makes 
it very similar to the template SN 1992A. Like many such 
events, SN 2004eo showed a slow evolution of the Si n line ve- 
locity l|Benetti et al.ll2005l ) , measurements of which indicate 
that Si extended down to layers at ~ 8500 kms" 1 , which is 
deeper than in SN 2002bo. The earliest spectra show a mild 
HVF in the Ca 11 IR triplet, and possibly in Si 11 6355 A. 

The s pectra and the light cur ve of SN 2004eo were pre- 
sented by IPastorello et ail l|2007i). The sample, although 
temporally not as dense as for SN 2002bo, covers both the 
peak and the nebular phases, and is probably the best so 
far for a SNIa of this luminosity group, offering a complete 
view of the ejecta of SN 2004eo except perhaps for the high- 
est velocities. 

In Section 2 our modelling is presented, and its 
limitations assessed. In Section 3 the modelling of the 
photospheric-epoch spectra is described, while the modelling 
of the nebular spectrum is discussed in Section 4. In Section 
5 the abundance stratification obtained from the models is 
presented and discussed. In Section 6 we present and dis- 
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cuss a theoretical light curve computed on the basis of the 
abundances that have been derived. Finally, in Section 7 we 
give our conclusions. 



2 METHOD 

Our modelling strategy of is to follow the evolution of the 
SN in the early phases with a photosphere-plus-expanding- 
envelope code based on the Sobolev approximation, and to 
explore the innermost part of the ejecta modelling a sin- 
gle nebular spectrum taken ~ 1 yr later with a non-local 
thermodynamic equilibrium (non-LTE) nebular code. 

For the early-time spectra, we used a Mont e carlo code 
based on th at first descr ibed bvlMazzali fc Lucvl (|l993l ) and 
updated by iLucvl (|l999l l and iMazzalil (|2000h . The code as- 
sumes a sharp photosphere, where the luminosity is emit- 
ted, and adopts the Sobolev approximation to follow energy 
packets in their propagation through the SN ejecta, which 
are represented by a density-velocity structure. The emerg- 
ing luminosity L, the velocity of the pseudo-photosphere v p h, 
and an epoch t since the explosion are required input, as is 
a set of depth-dependent abundances. Energy packets rep- 
resenting photon batches can interact with the ejecta via 
either line absorption or electron scattering. Upon absorp- 
tion, a packet can be reprocessed into a different line, thus 
describing the process of line blocki ng which is es sential for 
spectrum formation in Type I SNe (jMa zzali 2000). No con- 
tinuum processes are assumed to occur above the pseudo- 
photosphere. The properties of the radiation field and of the 
gas are iterated to convergence, temperature structure is 
established in radiative equilibrium and the ionization and 
excitation structures are computed using a modified nebular 
approximation that describes deviat ions from local thermo - 
dynamic equilibrium (for details, see lMazzali fc Lucy| [l993). 
We use the code version that employs abundance stratifica- 
tion, as discussed in Paper I. 

The nebular spectrum has been modelled using a code 
that follows in ID the propagation and the deposition of the 



7-rays and positrons produced in the decay of 56 Ni into 56 Co 
and h ence 56 Fe using a Montecarlo scheme (|Cappellaro et all 
1997). Gas heating by the fast particles created upon y- 
ray and positron deposition and the cons equent coolin; 
via line omission arc described in non- LTE (lAxelrodlll98i 
iRuiz-Lapuente fc Lucvl 1 19921 ; iMazzali et ail 120011). Stratin- 
cation in density and abundance is used, as in lMazzali et al.l 
|2007br i. 

Before we describe our results in detail, it is necessary 
to discuss the possible limitations of this method. Extract- 
ing the element distribution from a series of SN spectra 
is an ill-posed inverse problem: too many of the parame- 
ters that are used to compute the synthetic spectra are not 
known. In particular, assuming a certain density distribu- 
tion artificially narrows down the number of possible solu- 
tions. Here we adopt the density-velocity distribution of W7 
(|Nomoto et al.l fl984). for both the early and late-time mod- 
elling. Although W7 describes a more luminous SN la than 
SN 2004eo, it represents a typical density structure. Possi- 
ble uncertainties related to using that particular model are 
discussed in Section 7. On the other hand, using any other 
explosion model is no better method, as the uncertainties re- 



lated to the particular model would simply compound with 
those related to the spectrum synthesis codes. 

Additional caveats stem from the design of the code that 
we use to model the photospheric-phase spectra. The use of 
a sharp photosphere is very effective, and it keeps the anal- 
ysis free from the details of the explosion model, but it does 
not correspond to reality, in particular at times after max- 
imum, when the photosphere resides within the 56 Ni zone 
and consequently a significant fraction of the 7-ray deposi- 
tion occurs above the photosphere. These photons may alter 
both the flux level and the line profiles. Fortunately, this is 
mitigated by the decreasing density, which guarantees that 
even at advanced stages most deposition still occurs below 
the photosphere. 

Additionally, since the nebular spectra do not suf- 
fer from this problem, and since the velocities sampled 
by the nebular lines exceed 7000kms _1 , combining early 
photospheric-epoch results for the outer ejecta and late neb- 
ular ones for the inner ejecta we effectively avoid being too 
affected by the lesser accuracy of the post-maximum models. 

Another shortcoming is the fact that we use a one- 
dimensional approach, while in reality SNe la ejecta may 
be affected to some extent by aspherical effects. However, 
three-dimensional spectral modelling can only be attempted 
starting from an explosion model, or the number of free pa- 
rameters would be impossibly large, and it is a path that 
runs contrary to our approach. 

When modelling spectra, it is practically impossible to 
establish objective criteria with which to judge the good- 
ness of a fit. In fact, it is not fair simply to measure the de- 
viation from the observed flux, because not all wavelengths 
have the same weight when it comes to deciding whether the 
physical quantities that determine the spectrum are prop- 
erly described. Deviations in wavelength are if possible of 
greater importance, because they signal an incorrect abun- 
dance, ionization or density distribution. Deviations in flux 
can be due to an incorrect luminosity or temperature, but 
also to incorrect line strengths, and these are again deter- 
mined by density, ionization, and abundances. Therefore, a 
human observer is probably the best judge of the goodness 
of a fit. 

In the case of a tomography experiment, it is not enough 
to fit a single spectrum, but the entire series must be repro- 
duced by a single set of input data. There may be solu- 
tions that fit a single spectrum better than the one finally 
adopted, but the overall result may not be as good as that 
obtained with a different set. Therefore, a double compro- 
mise must be reached, on a set of input parameters that give 
both acceptable individual fits and on a fair reproduction of 
the spectral evolution. 

Further details about the accuracy of our results are 
given in Sections 5 and 6. 



3 PHOTOSPHERIC-EPOCH SPECTRA 

We have modelled 8 early-time spectra of SN 2004eo, cover- 
ing epochs from —11 to +2 2 days relative to B m aximum. 
There are later spectra (see iPastorello et al.ll2007h , but we 
did not model them because the Montecarlo code becomes 
progressively less accurate at advanced epochs and because 
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they probe very low velocities that are covered by the neb- 
ular spectrum. 

We adopted a bolometric rise time of 19 days for 
SN 2004eo . This is somewh at less than the average £?-band 
rise time l|Riess et al.lll999h . as motivated by the fact that 
dimmer SNela seem to rise somewhat more rapidly than 
brighter ones. This is the consequence of the smaller opac- 
ity and the shorter photon diffusio n time in the less nuc learly 
processed ejecta of dimmer SNe ijMazzali et al.l lioOl'l. The 
main input values for the synthetic spectra are recapped in 
Table 1, and the synthetic spectra are shown and compared 
to the observed ones in Figures 1 and 2. Using this approach, 
we map the abundance distribution in the ejecta. First we 
briefly discuss each spectrum in turn. 



3.1 19 September 2004: -11 days 

The earliest spectrum of SN 2004eo has an epoch of ~ 8 
days after explosion (Figure la). Although the SN was not 
a very luminous one, at this early epoch the small size of 
the ejecta results in a relatively high temperature, as shown 
by the fact that the two Fe features in the blue, near 4300 
and 4800 A, are dominated by Fe ill lines. Other strong fea- 
tures are due to Ca II, Si II, S II, O I. The ratio of the two 
Si II lines (A5972 and 6355 A), which is a good temperature 

m 



indicator ([Nugent et al.l ll995) is rather large. As shown by 
Hachingcr et al. I (|2008h . this is not directly the effect of a 
low temperature, as then the ratio would show the oppo- 
site trend, nor i s it c aused by Ti 1 1 lines, as claimed by 
iGarnavich et all (|2004 ). as in that case the much stronger 
Ti II lines in the blue would completely change the spec- 
trum. It is the result of saturation of the 6355 A line as Si II 
becomes more abundant with respect to Si ill going to lower 
luminosities, hence lower temperatures. So this is only indi- 
rectly the effect of temperature. In the UV, the spectrum is 
blocked by Fe II, Fe III, Ni II, Ti II, Ti III, Cr II, Cr II and 
Mg II lines. The flux redistribution caused by line blocking 
determines the properties of the optical spectrum. 



3.2 24 September 2004: -6 days 

The second spectrum has an epoch of 12.7 days after ex- 
plosion (Figure lb). The increase in luminosity does not 
keep up with the increase in size, and the temperature is 
slightly lower. Most lines and their strength ratios are how- 
ever similar to the previous epoch, since the temperature 
drop is small. The most noticeable change is the increased 
strength of the S II complex near 5500 A (AA5433, 5454 A for 
the bluer feature, AA5606, 5640 A for the redder one). This is 
indicative both of recombination and of an increased sulphur 
abundance as deeper layers are exposed. The abundance of 
silicon is also higher near the photosphere, but this has less 
of an effect on the Si II lines, which were already quite strong 
in the earlier spectrum. In both this and the previous spec- 
trum, the strength of the O 1 7774 A line indicates that a 
significant fraction of the outer layers was left unburned or 
only partially burne d, as indicated by the absence of carbon 
lines (|Mazzalill20"0lh . 



3.3 27 September 2004: -3 days 

In the third spectrum, the last before B maximum with 
an epoch of 16.4 days after explosion, the trend towards 
lower temperatures continues, while again the most appar- 
ent evolution is in the strength of the S II lines (Figure lc). 
The IME abundances are not much higher at this depth 
(9500 kms -1 ), and reach a maximum here. The features 
that characterise the spectrum are the same as at the previ- 
ous epoch. In all the pre-maximum spectra the near-IR con- 
tinuum is reproduced reasonably well by the model, which 
never overestimates the flux by more than 20%. This indi- 
cates that the opacity, which is mostly due to lines, is suf- 
ficiently large at optical wavelengths that the assumption 
of an underlying black body is not grossly incorrect. Line 
blocking in the UV is stronger in this spectrum, as the pho- 
tosphere moves inwards and more metals can absorb UV 
photons. 



3.4 2 October 2004: +2 days 

The spectra soon after maximum show subtle but significant 
differences with respect to the pre-maximum ones. In partic- 
ular, starting with the spectrum of 2 October 2004, which 
has an an epoch of 21.4 days after explosion (Figure Id), 
Fe II dominates over Fe III. This is clearly seen in the shape 
of the feature near 4800 A, which now shows strong Fe II 
multiplet 48 lines (4923, 5018, 5169 A) and has the charac- 
teristic shape of spectroscopically normal but relatively dim 
SNe la. Also, the O I 7774 A line became weaker, indicating 
that the oxygen content near this rather deep (9000 kms -1 ) 
photosphere is small or zero. The ionization of other ions 
that determine the UV also changes, as it does for iron, so 
that now Ca II, Ti II and Cr II, together with Fe II, determine 
the UV opacity. On the other hand, this spectrum still has 
a rather blue optical pseudo-continuum. 

The model reproduces the observed spectrum at least 
bluewards of ~ 6500 A, with a near-photospheric composi- 
tion that is still dominated by IME, but where 56 Ni is in- 
creasing. The mass enclosed below the photosphere is now 
~ O.8M0, and we are beginning to see the bulk of the NSE 
material. The model now overestimates the continuum in the 
red, indicating that the assumption of an underlying black 
body is less appropriate at this epoch. The observed pseudo- 
continuum redwards of the Si 11 6355 A line must be caused 
mostly by flux redistribution to the red from line processes 
that take place in the blue and UV. Thus, the luminosity 
used for the model is probably higher than the real value. 
Nevertheless, the strength of the synthetic O 1 7774 A line is 
comparable to that of the observed feature. This is impor- 
tant in order to estimate the oxygen abundance. In general 
the abundances are not very different from those of the last 
pre-maximum spectrum. 



3.5 7 October 2004: +7 days 

The next spectrum (Figure 2a) has an epoch of 26.4 days 
after explosion, and it continues to show the evolution to- 
wards lower temperatures (a larger ratio of the two Si 11 
lines, the growth of the Si 11-Fe 11 feature near 4800 A) that 
started from the earliest phases. The iron absorption is now 
much stronger, indicating that we are beginning to see the 
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Table 1. Parameters of the synthetic spectra. 



date 


epoch 


logL 




T s 




(days) 


[erg s" 1 ] 


(kms- 1 ) 


(K) 


19 Sept 2004 


8.0 


42.66 


12000 


13791 


24 Sept 2004 


12.7 


43.03 


10400 


13926 


27 Sept 2004 


16.4 


43.07 


9500 


12424 


2 Oct 2004 


21.4 


43.10 


9000 


10621 


7 Oct 2004 


26.4 


43.00 


7300 


10060 


11 Oct 2004 


30.4 


42.85 


6000 


9499 


14 Oct 2004 


33.4 


42.74 


4700 


9746 


21 Oct 2004 


40.4 


42.53 


2800 


10367 




Figure 1. Fits to the earliest 4 spectra of SN 2004eo. Basic input parameters are shown and line identifications are marked. The spectra 
are as follows: a) (top left): t = —11 days; b) (top right): t = —6 days; c) (bottom left): t = —3 days; d) (bottom right): t = +2 days. 
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4000 6000 8000 4000 6000 8000 

rest wavelength [A] rest wavelength [A] 

Figure 2. Fits to the spectra of SN 2004eo in the post-maximum phase. Basic input parameters are shown and and line identifications 
are marked. The spectra are as follows: a) (top left): t = +7 days; b) (top right): t = +11 days; c) (bottom left): t = +14 days; d) 
(bottom right): t = +21 days. 



inner, fully burned zone. Placing the boundary of the tran- 
sition between incomplete (IME-dominated) and complete 
(NSE-dominated) burning zones somewhere between 7000 
and 9000 kms -1 one derives, using W7, an upper limit to 
the NSE mass of ~ 0.7 Mp>. Th is is in agreement with the 
results of iMazzali et all l|2007al ). 



3.6 11 October 2004: +11 days 

In this spectrum (Figure 2b), at an epoch of 30.4 days after 
explosion, the low temperature begins to affect the pseudo- 
continuum significantly: the shape of the spectrum is now 
much redder. Fe II lines have grown in strength. Not only 
is the Ni abundance higher, but 56 Co is decaying to 56 Fe. 
The fraction of stable Fe is also larger. This is a feature of 
SN 2004eo which is also borne out by the analysis of the 
nebular spectra (see Sect. 3). At this and later epochs, the 
feature near 5900 A is probably contaminated by Na ID. 



3.7 14 October 2004: +14 days 



This spectrum (Figure 2c) has a nominal epoch of 33.4 days 
after explosion, and is very similar to the one obtained three 
days earlier. Most features are the same as in the previous 
epoch. The Na iD line is growing in strength, as is typical 
especially of dim SNela. This feature is not reproduced be- 
cause of_the_wehJmownuncertainty in the sodium ionization 
(e.g. IMazzali et al.| [l997). Although the photosphere is now 
very deep (4700 kms -1 ), the Ca n IR triplet absorption 
seems to reach even lower velocities, possibly suggesting the 
need for a higher calcium abundance at low velocity, or an 
incorrect estimate of the ionization there. The composition 
is dominated by 56 Ni (60% by mass) and stable Fe (35% by 
mass), the rest being mostly silicon and sulphur. This agrees 
with the nebular spectrum, which samples the region of the 
ejecta that is near the photosphere at this epoch. 
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3.8 21 October 2004: +21 days 

The last of the early-time spectra we modelled has an epoch 
of 40.4 days after explosion (Figure 2d). It is very red, and it 
shows deep and broad absorptions. The strong emissions are 
still P-Cygni components caused by the absorption blends: 
the epoch is still too early for nebular emission lines to 
emerge, as indicated by the good fit to the Ca II IR emis- 
sion component. The broad features near 4500 and 5000 A 
are caused by Fe II lines, although some Ti II contributes 
to the bluer feature, as it does in very underluminous SNe 
l|Filippenko et al.lll992T ). The Na I D line is becoming very 
strong, while other lines are the same as at earlier epochs. 

The synthetic spectrum matches the observed one sur- 
prisingly well, considering the advanced epoch and the very 
deep photosphere (2800 kms -1 ), which suggests that the as- 
sumption of a grey photosphere is not unreasonable at this 
epoch. The reason for this is probably the fact that since 
iron dominates the composition near the photosphere, Fe II, 
Ti II, and, further to the red, Co II lines effectively block 
the radiation and form a pseudo-continuum throughout the 
optical domain. The composition obtained agrees with that 
derived from the late-time spectrum (Sect. 3). 



4 THE NEBULAR EPOCH 

One nebular-epoch spe ctrum of SN 2004eo was presented in 
iPastorello etall (|2007T l. It was obtained 227 days after B 
maximum, on 16 May 2006. We use this spectrum to inves- 
tigate the structure and composition of the inner ejecta. The 
m odelling result s were a lready included in the compilation 
of iMazzali et alj l|2007al ). but the spectrum has been recali- 
brated because of scatter in the photometry and the values 
are therefore slightly different. 

A one-zone version of the nebu lar code can be use d 
to determine basic parameters as in IMazzali et al. | (|199S| ). 
The intrinsic width of the emission lines required to match 
the broad blends is 7300 kms -1 , which is o n the low side 
of th e distribution of line widths in SNe la (Maz zali et al.l 
1998). Nevertheless there is good overlap in velocity space 
between the most advanced photospheric epoch spectra that 
we modelled and the nebular spectrum, ensuring complete 
coverage of the ejecta and the consistency of the analysis. 
Material at velocities above 8000 kms _1 makes a negligible 
contribution to the nebular spectrum, indicating that the 
densities are too low and that 56 Ni is not present in any 
large amount at those velocities. 

For the stratified modelling the density stratification of 
W7 is again adopted, and the depth-dependent abundances 
are modified in order to fit the spectrum. 

The abundances in the nebula are dominated by NSE 
material. The best fit (Figure 3) includes a 56 Ni mass of 
0.38Mq, located mostly between velocities of 2000 and 
8000 kms -1 . The abundance of 56 Ni reaches ~ 0.65 by mass 
between 3000 and 7000 kms -1 . The innermost zone is dom- 
inated by stable Fe and Ni. This is necessary in order to 
avoid the formation of sharp emission peaks for the strongest 
[Fe ill] and [Fe n] lines and to achieve an ionization balance 
that gives the correct ratio of the [Fe m]-dominated 4700 A 
line and the [Fe ll]-dominated 5200 A line. The mass of sta- 
ble Fe is O.27M , so that the total NSE mass is 0.65M o (the 
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Figure 3. The nebular-epoch spectrum of SN2004eo and a syn- 
thtic match. 



mass of stable Ni is very small, in order to avoid a strong 
[Ni n] emission near 7000 A. 

The strongest observed emissions are fitted quite accu- 
rately. Many of the weaker lines that are not reproduced are 
weak Fe or Co lines, for which the collision strengths are not 
well known. The poor reproduction of those features means 
that we may be underestimating the 56 Ni mass by ~ 5%. 
Apart from the strong [Fe n] and [Fe ill] lines, Na iD is 
clearly visible, as are Ca II lines, which are however strongly 
affected by [Fe n] emission. Mg i] 4700 A contributes to the 
blue wing of the 4700 A feature: a Mg mass of O.O4M0 is 
obtained from the fit. Silicon does not have strong lines in 
the optical domain, but it has a cooling effect since it pro- 
duces lines in the IR, notably at 1.6 microns. A low upper 
limit for the sulphur abundance can be established from the 
line near 4000 A. The oxygen abundance must be very low 
at velocities below 8000 kms -1 , or a strong [O i] 6300 A line 
would form. The carbon abundance must also be very low, 
or [C i] emission would affect the blue part of the Ca II IR 
triplet emission. The abundances above 8000 kms -1 are the 
same as those derived from the early-time modelling. Be- 
tween 2800 and 7000 kms -1 some difference exists, but this 
is never more than 10%. 



5 ABUNDANCE STRATIFICATION 

Connecting the results of early and late-time spectral mod- 
elling, we can reconstruct the abundance distribution as a 
function of depth in the ejecta. The results are shown in 
Figures 4a and 4b. The 11 most abundant species which 
contribute to the spectra are displayed. Minor species can- 
not be treated accurately owing to the lack of significant 
spectral features. Figures 4a and 4b show that oxygen dom- 
inates the outer ejecta, down to ~ 12000 kms -1 . This is not 
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a particularly deep oxygen zone, despite the overall dimness 
of SN 2004eo, and it is comparable to most other SNela 
ijMazzali et alJl2007al ). The absence of carbon is noticeable: 
no carbon features are present that can be modelled, and the 
upper limit as constrained by the early spectra is very low 



(a few percent in the outer zones). This is also a feature of 
most SNela, which may indicate partial burning of carbon 
to oxygen in the entire progenitor. 

IME dominate at velocities between 7000 and 
12000 kms -1 . Silicon is the most abundant IME, followed 
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Figure 5. a) (top) Sensitivity of the results to velocity. The red 
line shows the observed t = 26 A day spectrum and the blue line 
is the best fit. The black band shows spectra computed changing 
the velocity boundaries of the various shells randomly by up to 
f5% with respect to the best-fit values, b) (bottom) Sensitivity 
to abundances within fixed velocity shells: same as above, but 
changing the abundances by up to 25% within each shell. 



by sulphur and magnesium. The abundances of sodium and 
calcium are much smaller. IME are present out to the high- 
est velocities surveyed by our spectral models, at a level 
of ~ 10% by mass. They do not extend significantly below 
~ 6000 kms -1 , indicating again that the innermost part of 
the ejecta was completely burned to NSE. Among the NSE 
isotopes, 56 Ni is the most abundant species between 3000 
and 7000 kms -1 , extending out to ~ 12000 kms -1 , while 
stable NSE isotopes (mostly stable Fe) dominate the inner 
3000 kms -1 , as discussed in the previous section. 

We can compare this distribution with those of W7 
(Figure 4c) and of SN2002bo (Paper I). In SN2004eo NSE 
material does not extend as far out as in SN 2002bo, where 
the bulk reached ~ 10,000 kms -1 but NSE material was 
present out to 15,000 kms -1 at the 10% level, while W7 is 
somewhere in between and shows no mixing-out. The outer 
extent of the IME is similar in SNe 2004eo and 2002bo, con- 
firming the result of lMazzali et alj (|2007al ). A peculiarity of 
SN 2004eo is the relatively large stable NSE material inner 
zone. 

The abundance distribution we have presented is ob- 
tained from our best-fit models, as discussed in Section 2 



8000 



above. It is instructive to verify the effect on the spectral fits 
of adopting different parameters. In Figures 5a and 5b we 
show the variation induced in the synthetic spectra by ran- 
domly varying the velocity boundaries of the various abun- 
dance shells by up to 15% (Fig. 5a) and the abundances 
by up to 25% (Fig. 5b). A relatively small change in the 
position of the shells can affect the result significantly. The 
sensitivity to abundance within a fixed shell is smaller, and 
the uncertainty can be estimated to be ~ 25% in the photo- 
spheric epoch. The examples shown here address the issue of 
fitting a particular spectrum. When trying to fit the entire 
time series, many of the solutions that may seem accept- 
able based on a single epoch do not result in a consistent fit 
of the spectral evolution. The masses of the NSE elements 
are estimated mostly from the nebular model, and have a 
smaller uncertainty. 



6 A LIGHT CURVE MODEL 

In Paper I we could show that using the abundances de- 
rived from the spectral analysis to estimate the opacity in 
the ejecta a more accurate synthetic light curve could be 
obtained than any based on theoretical explosion models. In 
particular, the mixing out of 56 Ni is not predicted in most 
models, which therefore fail to reproduce the early rise of 
the light curve. Unfortunately, the situation is not quite as 
good for SN 2004eo, possibly because there are fewer very 
early spectra, which makes it difficult to constrain the abun- 
dances above 15000 kms" 1 with high accur acy. In Figure 5 
the b olometric light curve of SN2004eo l|Pastorello et all 
120071 ) is compared to a synthetic bolometric light curve ob- 
tained with a Monteca rlo code based on that described in 
ICappellaro et afl (|l997l ). In the calculation, the optical pho- 
tons produced by the deposition of 7-rays and positrons dif- 
fuse in the ejecta, where they encounter an optical opacity 
which is determined by the abundance of the different iso- 
topic species as follows: 



K opt = 0.1[M(NSE)/M Q ] + O.O1[M(IME)/M ] 



(1) 



The synthetic light curve based on the abundance distri- 
bution derived from spectrum synthesis (Fig. 5, solid curve) 
resembles the observed one, with a peak at ~ 19 days, but 
it lies below the "observed" one by ~ 0.2 mag at all epochs 
before ~ 30 days, although the tail is well reproduced. 

The error introduced by the assumption of a sharp pho- 
tosphere in the models for the photospheric-epoch spectra 
can be assessed by verifying in the synthetic light curve what 
fraction of the 7-ray deposition occurs above the momentary 
photosphere. For example, at maximum light less than 10% 
of the deposition occurs above the photosphere, which is 
then located at ~ 9000 kms -1 , above the bulk of 56 Ni. This 
fraction increases to ~ 50% by day 40, when the photosphere 
is well within the 56 Ni zone. However, the nebular spectrum 
covers velocities out to ~ 9000 kms" 1 , and so the uncer- 
tainty caused by the use of a sharp photosphere is limited. 

One of the possible reasons for the discrepancy betwen 
the synthetic and observed light curves is the estimate of 
the mass and distribution of 56 Ni. 

The synthetic light curve can be improved by chang- 
ing the 56 Ni mass and distribution. Some possibilities are 
shown in Figure 6. First, if the amount of 56 Ni in the outer 
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Figur e 6. The b olometric light curve of SN 2004eo 
llPastorello et al.|[2007l) compared to a synthet ic light curve ob- 
taine d using our Montecarlo light curve model jCappellaro et all 
Il997l) and the abundances obtained from the spectral fitting 
discused in this paper. Although the synthetic light curve 
resembles the data, it fails to reach the same peak luminosity. 
The inset is an enlargement of the epoch of maximum. 

shells was underestimated because of poor spectral coverage 
or uncertainties in the spectral modelling at early times, a 
different rising part could be obtained. Adding 56 Ni at ve- 
locities above 10000 kms -1 improves the rising part up to 
~ one week before maximum, but it still does not match the 
peak. In order further to improve the fit, it is necessary to 
play with the distribution of 56 Ni in a way that is not jus- 
tified by the spectroscopic analysis. A model with enhanced 
56 Ni in the region between 7000 and 10000 kms -1 (Fig. 6, 
dashed line) results in a better fit of the peak. This model 
has M( 56 Ni) = O.49M , and M(NSE) = O.77M . However, 
because of the very large 56 Ni mass this model is much too 
bright on the tail. 

One further epicycle is then to remove 56 Ni at low ve- 
locities, below 6000 kms -1 . This model (Fig. 6, dotted line), 
with a total 56 Ni mass of 0.45M Q and M(NSE) = O.76M , 
fits the entire light curve, but is not physically justified. For 
example, a nebular spectrum based on this distribution of 
56 Ni would not at all resemble the observations. 

The question is then what is the source of the discrep- 
ancy. There are at least 4 possibilities: 

1) Calibration of the nebular spectrum. Incorrect cal- 
ibration would result in an incorrect estimate of the 56 Ni 
content in the inner ejecta, and hence in an incorrect light 
curve. This is a possibility since there is scatter in the pho- 
tometry obtained near the epoch of the nebular spectrum 
l|Pastorello et al.l 120071 ). It is however unlikely to result in 
a different relative distribution of 56 Ni and stable Fe-group 
isotopes, as this would directly affect the ionization balance 
and hence the line strength ratios. 

2) Estimate of the amount of 56 Ni in the region between 
7000 and 10000 kms" 1 . The early-time spectra are relatively 
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Figur e 7. The b olometric light curve of SN 2004eo 
llPastorello et al.l 12007^ compared to two synthetic light curves 
obtained by arbitrarily modifying the mass and distribution of 
56 Ni. In the model shown as a dashed green line 56 Ni was added 
at velocities between 7000 and 10000 kms" 1 , for a total M( 56 Ni) 
= 0.49Mq. This model fits the peak of the light curve, but is 
too bright on the tail. In the model shown as a dotted red line 
the additional 56 Ni at intermdiate velocities was compensated 
for by removing 56 Ni at velocities below 6000 kms" 1 , for a total 
M( 56 Ni) = O.45M . This model provides a good fit to the light 
curve but it is not physically motivated. The model from the 
abundance stratification [M( 56 Ni) = 0.38A/q] is also shown for 
comparison (thin blue line). 



insensitive to variations of the Ni mass in this velocity 
range: these layers are visible above the photosphere only 
after maximum, when most 56 Ni has decayed to 56 Co but 
only a small fraction of this has decayed to 56 Fe. Fe lines 
are therefore not greatly affected, and the effect of replacing 
nickel with cobalt is small. Nebular spectra, on the other 
hand, are also not very sensitive to small variations of the 
Ni content at these velocities because of the low overall 
density. 

3) Estimate of the bolometric correction and its evolu- 
tion. Evidence for this may be the fact that even though the 
spectral fits at maximum require an input luminosity com- 
parable to or somewhat lower than the corresponding LboI, 
the synthetic spectra thus obtained match the optical part 
of the observed spectra but are too bright in the near IR. 
Thus, the real value of LboI at peak may actually be smaller 
than estimated in the bolometric light curve. 

4) Uncertain treatment of the opacity. This is not likely 
to be a problem at late epochs, when the ejecta are transpar- 
ent, but it may affect the light curve calculation near max- 
imum. On the other hand, our synthetic light curve tracks 
the observed one quite well, the only problem being that it 
lies always below the data. This suggests that if there is a 
problem with the opacities, it lies in the 7-ray opacity, but 
this is a quantity that is well known and is unlikely to be 
grossly in error. The part of the light curve where the opti- 
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cal opacity is actually very likely to be incorrectly treated 
is between roughly 50 and 100 days, when the SN is making 
the transition to the nebular phase. This however has no 
effect on the estimate of the mass of 56 Ni either at the peak 
or on the tail of the light curve. 

Alternatively, a direct estimate of t he mass of 56 Ni 
may b e attempted using the prescription of lStritzinger et all 
(2006): L = 2 x 10 43 M( 56 Ni). With a bolometric light curve 
peak of 1.15 x 10 43 ergs _1 , this yields M( 56 Ni) = O.57M . 
This seems a very large value, given the results of the vari- 
ous light curve tests and the nebular spectrum. If we include 
56 Ni in the inner part of the ejecta, as suggested by the spec- 
troscopic modelling and by most explosion models, we can 
accommodate at most 0.49Mq of 56 Ni. This yields a good 
fit to the peak of the light curve, but it is too bright on 
the tail and it cannot reproduce the nebular spectra. Ad- 
ditionally, in all light curve models we have presented, we 
obtain a relation L ~ 2.5 x 10 43 M( 56 Ni), which again sug- 
gests that M( 56 Ni) = O.45M . The higher efficiency of con- 
version of 7-ray energy to opti cal light in our models with 
respect to the assumptions of IStritzinger et all (|2006T ) de- 
pends essentially on adopting a realistic density profile and 
a distribution of 56 Ni in the ejecta. Given all of the above 
uncertainties, a conservative estimate of the 56 Ni mass of 
SN2004eo is M( 56 Ni) = 0.43 ± O.O5M . 



7 CONCLUSIONS 

Despite the relatively small number of spectra available for 
modelling, a thorough description of the composition of the 
ejecta of SN 2004eo was obtained. This SN had a small 56 Ni 
mass for a normal SN la. The exact value is made somewhat 
uncertain by difficulties with deriving the bolometric light 
curve as well as by uncertainties in the calibration of the 
nebular spectrum and the behaviour of the optical opac- 
ity. Still our estimate, M( 56 Ni) = 0.43 ± O.O5M , places 
SN 2004eo among the least luminous normal SNe la, along 
with the template SN 1992A. 

A possible caveat concerns the total mass of the ejecta 
and the energy of the explosion. The assumption of the 
Chandrasekhar mass and of the W7 model yields a rather 
sensible abundance pattern and reasonably good synthetic 
spectra, suggesting that any deviation from that model 
should not be very large. Had the mass been significantly 
smaller, in our W7-based models we would have had to 
"hide" mass in elements that make a small contribution to 
the spectra: most likely oxygen in the photospheric phase 
and silicon in the nebular phase. This was however not 
required, and most elements included in the mixture do 
contribute significantly to the spectra at the appropriate 
epochs, via absorption lines at early times or emission lines 
in the nebula r phase. As for the k inetic energy, applying 
the formula of IWooslev etall (|2007h . E K = [1.56M( 56 Ni) + 
1.74M(stableNSE) + 1.24MQME) - 0.46]10 51 erg, we ob- 
tain Ek = 1.1 ± 0.1 x 10 51 erg. This is somewhat less than 
the original W7 energy (1.3 x 10 51 erg), but not so much 
less: the smaller contribution of burning to 56 Ni is in fact 
compensated by the larger production of IME, which con- 
tribute to the kinetic energy only ~ 20% less per unit mass 
synthesised, and by the large production of stable NSE iso- 
topes, which actually contribute more to the kinetic energy 



than burning to Ni. The main uncertainty regarding the 
kinetic energy comes from the outer extent of burning to 
NSE and the degree to which the outer carbon and oxy- 
gen were burned to IME. The zone between 9000 and 11000 
kms~ 1 still contains some oxygen. A signature of the lack of 
incomplete burning products at high velocities may be seen 
in the lack of HVF's, at least at the epochs sampled by our 
spectra. 

In conclusion, within the uncertainties SN 2004eo was 
a Chandrasekhar-mass explosion of a CO white dwarf 
that produced a smaller amount of 56 Ni than average 
(~ O.43Af ). The amount of stable Fe-grou p material (~ 
O.25M ), is somewhat larger than average |Mazzali et al.l 
l2007al ). The outer extent of IME (~ 1.1M ) is similar to 
most SNe la [the e xception being those defined as High Ve- 
locity Gradient by iBenetti et al.1 (2005)], so that the IME 
production was larger than average (~ O.4M ). Since there 
is no indication of unburned carbon, most of the remaining 
~ QAMq of material should be unburned oxygen. 

Such a low luminosity explosion may be the result of a 
strong initial deflagration phase, during which NSE material 
is produced and the star pre-expands, so that if and when a 
delayed detonation occurs this is not very eff ective, burning 
mater ial at low density mostly to IME only l|Mazzali et al.l 
l2007al ). Uncertainties arise mostly from the spectral cali- 
bration, the bolometric correction and the lack of extremely 
early data, calling for more complete coverage of SNe la so 
that the parameters of the explosions can be derived with 
even greater precision. Nevertheless, SN 2004eo is an impor- 
tant piece in the description of the variety of SNe la and in 
the effort to understand their origin. 
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